Plasmon dispersion in metal nanoparticle chains from angle-resolved scattering 
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We present angle and frequency resolved optical extinction measurements to determine the disper- 
sion relation of plasmon modes on Ag and Au nanoparticle chains with pitches down to 75 nm. The 
large splitting between transverse and longitudinal modes and the band curvature are inconsistent 
with electrostatic near-field models. Good agreement with an electrodynamic point-dipole model 
implies lower propagation loss and a larger signal bandwidth and group velocity than previously 
predicted. We conclude that for the design of optical nanocircuits coherent far-field couplings across 
the entire circuit need to be taken into account, even at subwavelength feature sizes. 

PACS numbers: 42.25.Fx, 42.79.Gn, 78.67.Bf, 71.45.Gm, 73.22. 
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A fundamental limit to the realization of sub- 
wavelength (sub-A) optical devices is that the interac- 
tion strength of dielectric objects with light vanishes as 
the objects gets smaller [l|, |J. Plasmonics may allow to 
overcome this inherent limitation of dielectrics by pack- 
ing the large polarizability of free electron resonances 
into a small physical volume 0, 01 • In this framework 
plasmon particle arrays have been proposed as an ideal 
platform that combines the ease of controlled nanofab- 
rication with the prospects of creating, e.g., ultra-small 
antennas to efficiently harvest, enhance and emit opti- 
cal power as well as a toolkit for nanophotonic cir- 
cuits Thus plasmon chains may act as sub-A width 
waveguides, waveguide bends, signal splitters, and fre- 
quency filters 0, [1, 0, As a parallel development, 
sub-A arrays of scatterers with magnetic rather than elec- 
tric resonances have recently gained tremendous interest 
for developing optical metamaterials. 
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Pioneering experiments have focused on qualitative un- 
derstanding of the resonance splitting and mode struc- 
ture in plasmon particle clusters and arrays [3]. The ob- 
served polarization-dependent resonant splittings in lin- 
ear (ID) particle chains, for instance, correlated well 
with trends anticipated from a simple near-field quasi- 
electrostatic model for the chain dispersion relation [1, 
[lol |. This qualitative 'quasistatic' model, which is valid 
on deep sub-A length scales, also formed the basis for 
forecasting the functionality of more complex struc- 
tures, such as plasmon chain splitters and multiplex- 
ers S, [i3|. Very recently, however, several groups 
have developed electrodynamic models that predict large 
quantitative and qualitative deviations from the qua- 
sistatic insights [l3|, [3, [I^j [iBl- If these deviations in- 
deed occur, a redevelopment of general design rules for 
complex nanophotonic circuits is required to include elec- 
trodynamic effects even on small length scales. 

Sofar, quantitative experiments to discriminate be- 
tween quasistatic and electrodynamic predictions at sub- 
A spacings have not been reported. In this Letter we 
present angle and frequency resolved optical extinction 



experiments on many arrays of Ag and Au nanoparticles 
at various particle sizes and sub-A pitches down to A/5. 
Sub-A plasmon arrays are an interesting system to test 
for dynamic effects. Recently, large modifications were 
predicted for the dispersion of modes in ID plasmon par- 
ticle chains and, equivalently, in ID magnetic split ring 
resonator arrays [l7|. For instance, the deviations in- 
clude a much larger fc = splitting between the longitu- 
dinal and transverse modes, and a polariton splitting of 
the transverse dispersion branch at the cros sing with the 
vacuum dispersion relation (light line) [l4 , iBl ■ For 
sub-A waveguiding these models further imply a large in- 
crease in group velocity and decrease in propagation loss 
compared to electrostatic predictions [f6|. The disper- 
sion relations measured in our experiment agree well with 
the recent electrodynamic models, and deviate strongly 
from quasistatic predictions down to pitches as small as 
A/5. As a consequence we anticipate that design propos- 
als for sub-A nano-optical circuits at currently realistic 
sizes (~ 50 nm pitch) can not be based on quasistatic 
analysis, because coherences and coupling across the full 
structure will dominate the optical performance. 
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FIG. 1: (A, B) Scanning electron micrographs of particle ar- 
rays on glass with r = 25 nm, d = 75 nm (A, Au) and r = 25, 
d — 100 nm (B, Ag). (C) Overview of the angle-dependent 
transmission setup. A collimated beam illuminates a large 
sample area. The sample angle is varied. Light is collected 
from a small spot on the sample (collection NA= 0.2). 
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We determine the nanoparticle chain dispersion rela- 
tion above the hght hne by far-field extinction measure- 
ments on Ag and Au nanoparticle arrays prepared on 
glass using electron-beam lithography. After physical va- 
por deposition of Ag or Au and resist liftoff, we obtained 
linear arrays of particles of 50 nm height, at pitches of 
d — 75, 100, 120 and 150 nm and with radii varied be- 
tween r = 25 and 55 nm, as determined by scanning elec- 
tron microscopy (SEM, cf. Fig. [HA,B)). The estimated 
error in determining r is approximately 2 nm. For each 
pitch the particle radii are below r/d = 0.37. We have 
fabricated square fields containing many parallel parti- 
cle chains (chain length 60 /^m), in which the inter-chain 
spacing is varied randomly (minimum 700 nm, mean 1 
/um) . The large inter-chain spacing ensures that coupling 
between adjacent chains is negligible, while the random 
variation suppresses grating effects that occur in exper- 
iments for periodic arrangements. Finally, we spin-coat 
the samples with a 100 nm layer of PMMA to ensure 
that the particle chains are embedded in a homogeneous 
dielectric environment. To determine the dispersion re- 
lation we use a wavelength-resolved transmission setup 
(see Fig. [IJC)) in which the sample is mounted on a 
rotation stage (axis perpendicular to the chains, along 
the PMMA/glass interface) that allows access to all in- 
cident angles from 6 = —60° to -1-60°. The chains are 
illuminated by a well coUimated white-light beam (di- 
vergence ~ 5°) from a fiber-coupled incandescent source, 
which illuminates a large (mm-size) area on the sam- 
ple. Using a pinhole on the transmission side, only the 
transmitted intensity from a ~ 15 /im spot (associated 
Ak/k ~ 0.09 (i.e., 5°)) on the sample is collected by a 
cooled Si-CCD coupled spectrometer. The transmission 
is obtained by normalizing the transmitted intensity at 
each angle to that recorded from an unpatterned sub- 
strate at the same angle. Using a broadband polarizer, 
we select the polarization of the input beam to be either 
transverse to the chains, or longitudinal (p-polarization). 
As the sample is rotated, the p-polarization also acquires 
a field-component transverse to the chain. 

Figure [2] shows typical transmission spectra for the full 
angular range for Ag particle chains of pitch 150 nm, and 
particle radius 45 nm (Fig [21^ A, C)), and for a smaller 
pitch of 100 nm and radius 30 nm (Fig[2)jB,D)). At nor- 
mal incidence {9 — 0°), a band of extinction around 
A = 500 nm for d = 150 nm (470 nm for d = 100 nm) is 
observed for transverse polarization, and at A = 625 nm 
(550 nm for d = 100 nm) for longitudinal polarization. 
The redshift of modes for the larger particles is consis- 
tent with the well-known single-particle resonance shift 
with increasing particle size [3, ll8|. The occurrence of 
two bands shifted to either side of the single particle res- 
onance is consistent with earlier reports on the 9 = 0° 
extinction of nanoparticle chains [9|. Qualitatively this 
splitting corresponds to the excitation of collective modes 
in the chain of dipole scatterers: the transverse mode is 
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FIG. 2: (Color online) Transmission as a function of angle and 
wavelength for longitudinal (A,B) and transverse polarization 
(C,D) for Ag arrays with r — 45, d — 150 nm (A,C) and r = 
30, d — 100 nm (B,D). The color scale runs from Tmin = 0.65 
in (A,C) and T^in = 0.80 in (B,D). 



blue shifted due to the antiparallel orientation of each 
dipole with the field of its neighbors, while the longitu- 
dinal mode is lowered in energy due to the alignment of 
each dipole with the field of adjacent dipoles. As the 
angle of incidence is increased, the two branches display 
opposite curvature, both shifting towards the single par- 
ticle resonance. 

The spectral dependencies of the extinction branches 
evident in Fig. [2 are qualitatively consistent with both 
the quasistatic model and full dynamic calculations for 
the dispersion relation of plasmon chain excitations. In 
order to quantitatively compare the data with the two 
models, we determine the transmission minima from a 
Gaussian fit to the transmission spectrum (plotted in the 
frequency domain) for each angle. In figure [3l we plot 
the resulting center frequencies taken from Fig. [2l^B,D) 
as a function of |k|| | = 27r/Asin(|6'|), i.e. the wave vector 
component of the incident beam along the chain. The 
measurements reach up to fc|| = 0.86ci;/c, or up to 60% of 
the light line in the medium embedding the particles. To 
appreciate the large width of the extinction resonances, 
the bandwidth at 1/e height of the fitted Gaussians is 
shaded in the diagram. The typical 1/e full width is 
approximately 2500 cm^^ (~ 60 nm). First, we compare 
to the generic quasistatic point-dipole prediction for the 
nanoparticle chain dispersion relation [^, [lo| : 
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Here loq is the single-particle resonance frequency, and 
= 2 for transverse, and kl = —4 for longitudinal 
modes jT9l|. We assume spherical particles with r/d taken 
from SEM data to obtain the quasistatic prediction in 
Fig- 131 (dotted curves) . There are two striking discrepan- 
cies between the data and the quasistatic model. First, 
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the splitting at normal incidence (fc|[ = 0) between the 
two branches is a factor two to three larger in the data 
than in the quasistatic model. The quasistatic model 
even falls outside the broad width of the extinction peaks. 
Secondly, the quasistatic model predicts that the trans- 
verse and longitudinal branch cross at kd = OAQtt inde- 
pendent of r/d at a frequency equal to the single-particle 
resonance. No sign of this crossing is observed in the 
data. It seems surprising that such discrepancies between 
data and the quasistatic model haven't been noted in ear- 
lier studies tQi] . These studies focused on fc = only, with- 
out investigating nonzero scattering angles. On the basis 
of fc = data only, one might assume that a larger split- 
ting is due to an error in r/d. For our data, this would 
imply an unlikely 40% error in estimating r/d. However, 
even if one would scale r/d to match the = splitting, 
the quasistatic model would still not be consistent with 
the full angle-dependent data set: the presence of the 
band crossing at kd — 0.467r would impose a much larger 
curvature of both bands than observed in our data. 

The fact that the quasistatic dipole model does not 
describe the data could be due to several approxima- 
tions: The quasistatic model ignores dynamic effects, 
multipole effects, the presence of the PMMA-air inter- 
face, and particle anisotropy. Based on Ref. we 
conclude that multipole effects do not resolve the dis- 
crepancy between the quasistatic model and our data: 
multipole effects do not alter the quasistatic dispersion 
relation for r/d < 0.35, a criterion satisfied in our ex- 
periment {r/d = 0.3 ± 0.02). For larger r/d multipole 
effects in fact reduce the quasistatic k = splitting, and 
shift the band crossing to even smaller wave vectors [13] , 
inconsistent with the absence of a crossing in our data. 
To study the effect of the dielectric interface, we have 
analyzed the quasistatic model for dipoles in front of a 
dielectric interface on the basis of image dipole theory [l[ . 
Within this model the image dipoles have a weaker dipole 
moment by a factor {cm — l)/(em + 1) (with Cm the em- 
bedding dielectric constant) and are located at a distance 
of 150 nm, i.e., further than the array pitch. As a con- 
sequence the estimated increase in splitting is less than 
15%, i.e., much less than the 2 to 3-fold enhancement ev- 
ident in the experiment. Regarding anisotropy, the parti- 
cles for which data is reported in Fig. [3] arc mildly oblate 
(height 50 nm and diameter 60 nm). At 6* = 0°, both 
polarizations used in the experiment are along equiva- 
lent axes of the single-particle polarizability tensor, so 
that an enhanced splitting is not due to excitation of 
distinct single-particle resonances. In addition, we have 
calculated the polarizability tensor (including dynamic 
depolarization shifts I,!!!). The resonance relative to 
that of a r = 30 nm spherical particle is only shifted 
by 4 nm, and the on-resonance polarizability along the 
long axes is in fact 10% smaller. Within the quasistatic 
model, particle anisotropy is hence only expected to re- 
duce the splitting, rather than explaining the observed 
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FIG. 3: (Color online) Symbols: frequency (in cm~l (left 
axis) or eV (right axis)) of minimum transmission versus 
wave vector for transverse and longitudinal polarization in 
Fig. [21;B,D), i.e. the d = 100 nm, r = 30 nm Ag particles. 
Shaded areas: 1/e bandwidth of the transmission minima. 
Straight lines: light lines in air and in the embedding medium. 
Thick (thin) curves: dynamic (static) prediction. 



enhancement. 

Having excluded that modifications to the quasistatic 
model due to multipole effects, interface corrections or 
particle anisotropy can explain the dispersion relation 
observed in the experiment, we now compare the data 
to a full electrodynamic point-dipole model [lij that in- 
cludes Ohmic damping, radiation damping, and depolar- 
ization shifts in the single-particle polarizability [18], as 
well as all terms in the dipole field. Within this model the 
dispersion relation has recently been calculated pertur- 
batively [l3. a nd self-consistently for finite [iH] and infi- 
nite arrays [l6|. Using a dielectric model for silver that 
is a modified Drude fit [2^ to tabulated data 21], and 
taking r and d from SEM observations, we plot the self- 
consistent infinite array model as solid curves in Fig. [31 
This model predicts a band splitting at fc = that is ap- 
proximately twice the splitting in the quasistatic model 
for the same geometrical parameters, as well as a band 
curvature for wave vectors away from fc = that is 
in much better agreement with the experimental data. 
Keeping in mind that the dynamic model has no ad- 
justable parameters, the reasonable correspondence with 
the data is a strong indication that retardation effects 
and far-field coupling are indeed determining the disper- 
sion relation for d = 100 nm silver nanoparticle chains. 
The remaining deviation from the dynamic model may 
yet again be due to the small particle anisotropy, the 
PMMA-air interface, or multipole effects [l3[ the influ- 
ence of which is much more difficult to estimate quantita- 
tively for electrodynamic rather than quasistatic models. 

Based on the data in Figs. [2l [3] we conclude that the 
large deviation of the splitting between the transverse 
and longitudinal branch at normal incidence compared 
to the quasistatic prediction is a good indicator for the 
relevance of far-field effects. Figure 2] shows the measured 
relative splitting Auj/lo for many combinations of parti- 
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FIG. 4: (Color online) Symbols: Relative splitting between 
transverse and longitudinal branch at = versus r/d for 
several pitches d as indicated for Ag and Au arrays. Black 
curve: quasistatic prediction for Ag arrays. Colored curves: 
dynamic predictions for d = 75 nm (Ag, dashed) and d = 
150 nm (Ag, dotted). Multipole effects set in beyond the 
dashed vertical ,10] . 



cle radius r and pitch d plotted against r/d, both for Ag 
and for Au particle arrays. In the quasistatic limit, the 
splitting is simply proportional to {r / dY (see Eq. ([IJ). 
In the dynamic model, the splitting depends on r and d 
separately as is clear from the dynamic predictions for 
the smallest [d — 75 nm) and largest [d = 150 nm) pitch 
used in our experiments (dashed/dotted in Fig. [4]). Fig- 
ure m demonstrates that for Ag particles the splitting is 
generally a factor ~ 2 larger than the quasistatic the- 
ory predicts, which is in agreement with the dynamic 
predictions. For gold particles, the difference is not as 
large, which may explain why previous studies did not 
resolve deviations from quasistatic theory |9|. Still, the 
observed splittings for Au systematically exceed the qua- 
sistatic prediction. We attribute this smaller effect for 
gold to the much lower albedo (< 20% for r = 25 nm 
Au, compared to ^ 80% for Ag) 0]. Since low-albedo 
particles radiate less strongly, far-field corrections to the 
dispersion relation will be less important. 

In conclusion, the experiment reported in this Letter 
shows that the dispersion relation for plasmon modes on 
metal chains of sub-A pitch is strongly modified by far- 
field retarded interactions, even for pitches as small as 
d = 75 nm or ~ ,2^5. This experiment thus confirms 
recent models 14 , 1^ 3] that overturn the conventional 
view of propagation along metal nanochains [i,[l^. The 
next challenge is to benefit from appealing consequences 
of these recent predictions that occur below the light line, 
i.e., in the range only accessible to near-field microscopy. 
In this regime propagation distances up to 5-10 /im at 
group velocities around 0.3c are feasible [l6| . These prop- 
agation distances far exceed initial quasistatic estimates, 
which were based on the damping rate of single parti- 
cles [§, ^ . Since a large part of single-particle loss is by 



radiation into the far-field, most loss can be suppressed 
by far-field destructive interference. Since far-field cou- 
pling is expected to be stronger in 2D clusters [l3|, this 
work implies that for designing nano-optoelectronics in 
general, an essentially quasistatic electrical circuit-design 
approach Q is only applicable to structures < 50 nm, 
i.e., a scale at which fabrication of several coupled com- 
ponents is extremely challenging. Designing coupled pho- 
tonic structures at the current fabrication limit will al- 
ways require a fully electrodynamic optimization of the 
coherent coupling between all building blocks: even at 
sub-A pitch far-field interference effects are key to opti- 
mize the functionality and reduce the loss. 
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